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Abstract 

We study the spin-polarization mechanism in the highly dense nuclear matter with the 
relativistic mean-field approach. In the relativistic Hartree-Fock framework we find that 
there are two kinds of spin-spin interaction channels, which are the axial-vector and tensor 
exchange ones. If each interaction is strong and different sign, the system loses the spherical 
symmetry and holds the spin-polarization in the high-density region. When the axial-vector 
interaction is negative enough, the system holds ferromagnetism. 
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1 Introduction 



Resent discovery of "Magnetar" which is a neutron star with super strong magnetic field, 
seems to revive a big question on the origin of the strong magnetic field. Since there is spread 
bulk hadronic matter beyond the nuclear density inside neutron stars, it should be interesting 
to consider its origin in the context of dynamics of hadronic matter; e.g., if the spin-polarization 
of baryons are realized in nuclear matter, ferromagnetism may occur in neutron stars. 

For quark matter, one of the authors (T.T.) has recently indicated a possibility of spin- 
polarization of quarks interacting with one-gluon-exchange (OGE) interaction [0; the Fock 
exchange interaction between quarks has a role to align spins, which is similar to the electron 
system Q . Using this result he suggested that super strong magnetic field expected in magnetars 
may be explained if they are quark stars. There he also found relativistic effects give rise to a 
new mechanism for ferromagnetism, which is never appeared in the nonrelativistic case. 

As for the normal nuclear matter Pandharipande et al. [Q] have reported no possibility 
of stable spin-polarization within the non-relativistic framework; magnetic susceptibility never 
changes its sign within densities relevant for neutron stars. 

On the other hand Niembro et al. ^ have suggested a possibility of spin-polarized 
nuclear matter using the relativistic Hartree-Fock (RHF) approach Q, though spontaneous 
spin-polarization occurs at too high density. The results in Refs. |5|, |6| suggest that the rel- 
ativistic framework may be more favorite for spin-polarization than the nonrelativistic one. 
Checking their framework we find some problems about the calculation. First, they implic- 
itly defined the spin-polarization of the system by using the eigenstates of the spin operator 
^z{= ^ 1) = 75707'^ in the rest frame of each particle, and all the particles take the same 
eigenvalues in the fully polarized state. This may be a direct analogue of the nonrelativistic fer- 
romagnetism. However, the spin operator cannot commute with Hamiltonian in the relativistic 
theories and thereby we must carefully treat two polarization degrees of freedom for baryons. 
Actually their choice is not a unique choice; the spin states of all the baryons do not necessarily 
need to be the eigenstates of the spin operator even in their rest frames. Instead, the spin 
orientation should depend on the momentum of each particle. Thus we must consider the spin 
configuration of the spin-polarized state in the phase space. 

Secondly, they kept spherical symmetry during the formulation though the spin-polarized 
system may break it due to the existence of a specific direction (we can assign it as the z axis 
without loss of generality) along the nonzero magnetization vector in the spin-polarized state. 
The interaction-energy density has generally the additional momentum dependence besides that 
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from the propagator in the relativistic theories due to the lower component of the Dirac spinor. 
Since there should be appeared a characteristic direction in the spin-polarized state, there is 
another scalar product of momentum and the unit vector along the z axis besides that of 
momenta in the interaction energy density; the expression of the interaction energy density is 
no longer rotation-invariant and retains only the rotational symmetry around the z axis. This 
may in turn suggest a possibility of axially symmetric deformation of the Fermi sphere. 

Thirdly, they used in Refs.[^, ^ the o", cj, p and 7r-meson exchange for the RHF calculation. 
In reality these interactions should be considered to be the in-medium one, not the bare one, 
and we do not have any enough information for the two-body interaction which should be used 
there, especially for the spin-spin interaction channels. 

Putting these remarks aside, we know there has been no systematic and sufficient discussions 
on this topic in the relativistic many-body approach Q, particularly in view of the relativistic 
effects. For example, the usual calculation, either in the non-relativistic or relativistic framework, 
has been done under the assumption of spherical symmetry for the mean fields. In this paper 
we reexamine the spin-polarization of nucleon matter within the RHF approach, focusing on 
the breaking down of spherical symmetry and importance of the relativistic effects. Since it has 
been shown in Refs.|Q, |5|, |6| that the Fock exchange interaction plays a key role in the context of 
the spin polarization, we must treat this matter within the RHF approach. Then we take some 
choices of the spin-vectors dependent on the momentum of each particle, and study relations 
between nuclear properties and these choices. 

In the next section we give our formalism, where we clarify what kinds of interactions are 
effective for the spin-polarization and figure out the role of the spherical symmetry breaking. 
We also classify the various spin configurations in the relativistic formulation, whereas there is a 
unique choice in the nonrelativistic theories. Results of numerical calculations are given in Sec. 
3. Sec. 4 is devoted to summary and concluding remarks. 

2 Formalism 

In this section we briefly explain our formulation to describe the spin-polarized system. There 
should appear a special direction along the spin-polarization; it is defined to be oriented to the 
positive z-direction. Such a system breaks spherical symmetry while the axial symmetry around 
the z-axis is preserved. 

In the RHF framework the interaction energy density in the isospin saturated system is 
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generally written as 

emi = ]^ j -^^-^^d'^k[Tr{iS{p)]Vs{p-k)Tr{iS{k)} 
+Tr{iS{p)^^}Vvip - k)Tr{iS{k)J^'} 
+Tr{iS{ph5}Vp{p - k)Tr{iS{kh5} 
+Tr{iS{ph5lt.}'DA{p - k)Tr{iS{k)j5j''} 

+Tr{iSip)a^,}VT{p - k)Tr{iS{k)a^'''} ] (1) 

for the one-boson exchange type interaction, assuming no derivative coupling. Here S{p) is 
the nucleon propagator with momentum p, Da {a = S, V, P, A, T) is the linear combination of 
meson-propagators with the nucleon-meson couplings. Note that these interaction terms do not 
necessarily appear in the original Lagrangian, and that some of them are given by other ones by 
way of the Fierz transformation. When using the a- (scalar), lo- (vector) and ??- (pseudo-scalar) 
meson exchanges, for example, 

^^(^^ = -^ + ^5^A.(g)-i5SA.(g) + ^5^A,(g) (2) 



Vp{q) = -^ + ^5^A.(g) + ^gSA^(g) + ^5^A,(5), (4) 



^ .2 A („\ ^ „2 A f„\ , 1 „2, 



^A(g) = -5^A,(g) - -5^A^((7) + -g^A^(g), (5) 
VT{q) = \gl/:^.{q)-\gl\{q), (6) 



where ga (a=o", u and rj) is the nucleon-meson coupling strength, and A^ is the meson-propagator, 

^a{q) = ^r' — 2 {a = a,uj). (7) 

The first constant term g^jrr?^ {gt/'^t) (^v) indicates the Hartree direct contribution, 

and other terms are the Fock exchange contributions. Introducing other mesons or vertex form- 
factors to couplings do not change the form of eq. (|^) . 

In the spin-polarized system, the expectation value of the spin operator has a nonvanishing 
value, and thereby at least the axial-vector (A) and tensor (T) exchange terms survive in eq. 
(|l]) as well as the scalar (S) and vector (T) terms. Then the nucleon self-energy must include 
their contributions, and the ground state breaks parity and spherical symmetries. 

In the relativistic mean-field (RMF) approach we usually neglect the momentum dependence 
of the propagator because the nucleon-nucleon interaction can be effectively treated as the zero- 
range one in low energy phenomena as far as the typical energy and/or momentum scale is much 



less than the meson mass. Actuahy only very small momentmn dependence has appeared in 
the full RHF calculation [Q, Here, we take the zero-range approximation for two nucleon 
interaction which can be described as follows: 

= ^ (« = '5,F,p,Ar). (8) 

Even with this approximation the RHF calculation is still complicated because there appear 
the axial-vector and tensor mean fields. These parity violation terms largely mix the positive- 
energy and negative-energy states in the single particle wave- function. Thus completely self- 
consistent calculations become very complicated and must be done very carefully as for the 
vacuum polarization. Furthermore we do not have sufficient information of all channels of 
the in-medium interaction between nucleons, especially in the axial- vector and tensor channels. 
Thus we will not get any clear conclusion in spite of very tough calculations by solving the RHF 
equation self-consistently. Instead of solving the exact self-consistent RHF equation, we take a 
variational approach in the RHF framework. 

The two degrees of freedom of the spin polarization for each nucleon is denoted by C = 1 and 
= — 1, which we call spin- up and -down, respectively. Then we take the nucleon propagator 
with four-momentum p in the following form; 

S[vX) = Sf{pX) + Sd{p.C) (9) 
with the propagators of a vacuum piece {Sp) and a density-dependent piece (S^), 



Sd{pX) = {l,P*^' + )} + |^n(p, 05{p, - e,), (11) 



where p*^' = p^' - U^^, M* = M - Us and E* = y/p"^ + M*2. In these equations Us and C/^ 
are the scalar and vector mean fields, and Sp is the single particle energy of the nucleon with 
momentum p, Ep = E* + Uq. a{p* , (") is the spin- vector of the nucleon with momentum p which 
satisfies the following conditions: 

apa" = -1, af.p*" = 0. (12) 

In the following we only keep the density-dependent piece S^) with the momentum distribution 
function n{p, C) to be determined, for which we assume the axial symmetry along the spin- 
polarization. Note that this form of the propagator is the same as the one when the nucleon 
self-energy includes only scalar Us and vector U^ mean fields, which are independent of the 
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nucleoli's spin. Moreover, we can easily see that the expectation value of the pseudoscalar 
operator automatically vanishes, 

Tr{iS{p)^^} = 0. (13) 

So the pseudoscalar (P) term in Eq. (|l]) vanishes in this case. 

Then the total energy density e-r is separated into two parts: the spin-independent part esiD 
and the spin-dependent part esn as 

= ^siD + <^SD- (14) 

The axial- vector and tensor exchange channels contribute to the spin-dependent part eso while 
the kinetic energy and contributions from the scalar and vector channels are involved in the 
spin- independent one esio- 

Under the zero-range approximation the spin-independent part of the Fock contribution 
can be incorporated into the Hartree one. Then it is possible to redefine the two-body scalar 
interaction by taking into account the Fock terms, which corresponds to the usual relativistic 
Hartree (RH) approximation. Thus we can write the scalar mean-field by the expectation value 
of the a field {a) as 

Us = gAa), (15) 

and calculate the effective mass M* with the usual RH approximation. The expectation value 
(cj) is given by the equation, 

d ~ 

d{(j) Y 
where the scalar density is defined as 

and U[a] is the self-energy potential of the sigma-field, whose expression is given in Ref. [0, 

In the RH approximation, furthermore, the spatial components of the vector mean-field are 
vanished, and the time component is given as 

Uo = gu^PB = 9lo^Pb{0, (19) 

where psiC) is the baryon density contributed from nucleon with the spin suffix ( as 

PB{0 = 2j (03^p;C). (20) 



Under this formulation the spin- independent part of the energy density eg id becomes 

^ . (27r)3 ^ 
while the spin-dependent energy density esD is calculated to be, 



esiD = 2Y: J ^nip;C)E; + U[a] + ^[J2pB{0r, (21) 



^SD = ^Pa + ^Pt (22) 
with the axial- vector and tensor densities, 

PA = I^Jr{SD{pHl'} = PB<^.>=Y.J^^(p-^O^a. (23) 

PT = J^,Tr{SD{p)<7i2} = PB<P^.>=J2j^n{p-C){a.-^^ao} (24) 

Other components of the axial-vector and tensor densities are vanished because of the axial 
symmetry of the momentum distribution. Note that the axial-vector and the tensor interactions, 
even though they are not necessarily included in the original Lagrangian, may arise from the 
Fock exchange interactions by way of the Fierz transformation, as is seen in Eqs. (5) and (6). 
In this sense we can say that the Fock exchange interaction is essential for the system to be 
ferromagnetic [Q, ^, ^ . 

In order to figure out the properties of the spin-polarized matter, we solve RH equation 
([l^ and calculate the energy-density by fixing the baryon density pB and the spin-polarization 
parameter Xg defined by 

Xs = (pT - Pi)/pB- (25) 

where pi = pb{C = ^) and pi = pb{C = ~1)- (For convenience the spin-up and spin-down states 
are indicated by the symbols j and J,, respectively.) 

In our approach the wave-function is not an exact solution of the Dirac equation in the 
mean-fields. So we need to specify the spin configuration in the system by choosing a spin- 
vector a^. Once it is fixed, the momentum distribution of the single-particle state should be 
also determined accordingly. The best way to this end is to choose the configuration to optimize 
the total energy density ey, while it may be rather complicated. Instead, we examine here the 
following three choices by physical considerations. 

The total spin-polarization is directed to the positive direction of the z-axis; here we define 
a unit vector Cm = (0,0, C)- Usually the spin- vector is chosen as {0,Cm) at the rest frame 
of the nucleon, we should call this choice as Choice-l(Chl). This choice may be a natural 
extension from the nonrelativistic ferromagnetism and has been also taken in the context of 
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ferromagnetism of quark matter Then the spin- vector with momentum p becomes 

_ ,^. (Cm • P)P 1 _ Cm • P /„„N 

"^-1^^+ M*{E; + M*y ' " - M* ' ^ 
by way of the Lorentz transformation. In this choice pA and pT can be written as 

PA = ^EC{2p.(C) + Pb(C)}, (27) 
C 

PT = ^EC{P.(C) + 2pB(C)}. (28) 

C 

Then the interaction-energy density can be written only in terms of the scalar Ps{C) and vector 
densities Pb{C) of nucleons with the spin C, and the expression still holds spherical symmetry. In 
the high-density limit, as mentioned before, the scalar density ps approaches to a finite value in 
the RMF approach |8|, and its contribution becomes negligible, so that we can see that p\ < p\- 
in the limit ps — > oo. Hence the spin-polarization occurs when Ca + 4Ct < 0. The above choice 
does not necessarily lead to the state with the spin-alignment, and the total spin per nucleon 
converges to one-third (< S^/A >= pa/ PB Xs/'i) in the infinite density limit; magnetization 
is not so large even if all the nucleon spins align. 

Here we consider another choice, Choice-2(Ch2) for the vector to give the maximum for 
\az\ within our framework. In this new choice the spin- vector becomes 

E;{Cm-p) 

"° " MV(Cm-p)^ + M*^' ^'^^ 

M^^Cm + (Cm • P)P 
MV(Cm • P)' + Af *^ 

Substituting the above form into eqs.(|23|) and (p^), we get 

PA 



2i:cy (^-(p,c)^^ — (31) 



PT 



^ f <fp , , M* 
2i:c/p^n(p;0^5j^ (32) 



When pb — > oo, px converges to a finite value, and pA is proportional to the baryon density 
Pb- If the momentum-distribution n{p;C) is taken as usual Fermi-distribution, the total spin 
per nucleon converges 

< EJA >= pa/pb^Xs/2. (33) 

Hence we can expect a ferromagnetic phase as long as Ca < 0. We shall see that this choice is 
the most appropriate for ferromagnetism. 
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As another choice, Choice-3(Ch3), we take the vector to give the maximum of the ex- 
pectation value < fiYjzjA >= pt/ PB-, which is reduced to the same value as that of the spin 
operator in the nonrelativistic limit. In this case the spin-vector becomes 

ao = C, (34) 



« = (35) 

with 

PT=P-CmiCm.-p)- (36) 

Then the axial-vector and tensor densities can be written as 

PA = 2Ec/7|t3-(p;C)^=^— (37) 



^ = ^Ec/^-teO^^^^. (38) 

In this choice pA and pT show the opposite behaviors in the infinite density limit; p^ converges 
to a finite value, and px is proportional to the baryon density pB- Hence the system becomes 
spin-polarized state when Ct < 0. However, the expectation value of the spin operator has 
a nonvanishing value and becomes 

< >= pa/pb^O, (39) 

which means that the system is not ferromagnetic in a usual sense. 

Since the spin-dependent energy esD depends on the coupling strengths, Ca and Ct, we have 
to carefully determine which choice is most appropriate. In the low density region pA px ~ 
XsPBi the spin-saturated system should be stable and then the two coupling constants Ca and 
Ct must satisfy the following relation 

Ca + Ct> 0. (40) 

If Ca > and Ct > 0, the spin-saturated system must be stable in all the density region, and 
spherical symmetry is always held. Otherwise the analysis at the infinite density limit show 
us that system becomes spin-polarized if Ca < < Ct (Ch2) or Ct < < Ca (Ch3) holds. 
Even in Chi the spin-polarization can occur if Ca + 4Ct < 0, but this condition expects large 
negative Ct and the spin- vector of a nucleon prefers Ch3 than Chi. 

Furthermore we should note that the above expressions of eq.([3l|) and eq.(p2|) do not preserve 
spherical symmetry. The relativistic effects automatically give rise to the spherical symmetry 
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breaking. From this fact we can naturally expect that the momentum distribution n{p;() is 
allowed to be distorted while keeping the axial-symmetry. In order to estimate the effects of 
distortion of the momentum distribution, we introduce the quadrupole-distorted distribution 
function n{p) as 

C) = no(e^(^Vx, e^^'^^py, e'^^^^V,; (), (41) 

where no{p; C) = ^(pf(C)~IpI) with the Fermi-momentum p^r. The parameter A(C) is determined 
to give the energy minimum of the spin-polarized system. 

As mentioned in the previous section, the two contributions from the axial- vector and tensor 
channels are not considered in the original interactions; instead they are derived by the Fierz 
transformation from the Fock exchange interactions in other channels. Using several parameters 
given in previous works, however, we get various spin-properties even at the normal nuclear den- 
sity. For example the parameter-sets of HF-I in Ref. within the a- and uj- meson exchanges, 
gives us Ca = —8.55 and Ct = 30.2, and the parameter-sets with vr-, a-, p- and 7r(PS)- meson 
exchanges in Ref. gives us Ca = 3210 and Ct = —3200. The latter extraordinary value 
comes from the pion-exchange {M'^g^/m^ = 8200), which never allow the zero-range interaction 
approximation. Thus the values of two coupling strengths, Ca and Ct, are still very ambiguous, 
and they cannot be individually determined at present both in theoretical and experimental 
ways. Instead of using parameters given in previous papers, we investigate the spin-polarization 
of nuclear matter by varying the values of Ca and Ct in this paper. 

As discussed above, we have a possibility to get ferromagnetic matter only in the case of 
Ca < < Ct, where Ch2 for the spin- vector must be most appropriate. Hence we study the 
spin-properties only in this case except in some figures where properties are also calculated with 
Chi for comparison. We shall see that ferromagnetism can occur due to the spherical symmetry 
breaking in the relativistic framework, by concrete numerical calculations in the next section. 



3 Results and Discussions 

In this section we make concrete calculations as for the spin-polarized nuclear matter and 
discuss their consequences. We calculate physical quantities concerning the magnetic properties 
only for the isospin symmetric matter and make comparison with neutron matter in the final 



place. In the actual procedure, first, we evaluate the sigma-field with eq.(16) by fixing baryon 
density ps and the spin-polarization parameter Xg = (/0| — Pi)/ Pb] namely we solve the RH 
equation. Secondly, we substitute the result into eqs. (pTI) and (|3^, and obtain the spin- 



dependent energy ego in eq. (|22D . Repeating these processes by varying A| = \{C, = 1) and 
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A| = A(C = — 1), we search the energy minimum. 

The parameter-sets PMl {M*/M = 0.7) and PM4 {M*/M = 0.55) are used for the RH 
calculation; the definite values of the parameters are given in Table 1. In Fig. 1 we show the 
density-dependence of the total energy per nucleon {Et/A) and the effective mass normalized 
by the bare mass (M* /M). The values of the coupling strengths in the axial-vector and tensor 
channels must be consistent with the spin-properties at normal nuclear density such as the 
spin-symmetry energy. We here define the spin-symmetry energy, inversely proportional to the 
magnetic susceptibility by 

Ssp-sym - ^2 l^==0- l^^J 

In this work we take its value as Esp-sym = 25(MeV), while it is not clearly determined from 
experimental information and we use three kinds of the parameter-set: Ca = (SDl), Ca = 
—50 (SD2) and Ca = —100 (SD3). Here we restrict ourselves to the cases with Ca < , since 
matter would be ferromagnetic only in this case. The detailed values of parameters are given in 
Table 2. 

In Fig. 2 we show the density-dependence of the spin-symmetry energies Egp-sym with PMl 
(a) and PM4 (b). If this value becomes negative, the spin-saturated system becomes unstable 
and the spin-polarized one is favored. The long-dashed, dashed and solid lines indicate results for 
SDl, SD2 and SD3, respectively. For all the parameter-sets the spin-symmetry energy increases 
in the low density region as baryon density becomes larger. While the spin-symmetry energy 
monotonously increases in the case of Ca = (SDl), it decreases and becomes minus above 
a critical density pc for the cases of negative Ca (SD2 and SD3): Pc/po = 8-74 for SD2 and 
Pel = 4.28 for SD3. Note that the spin-symmetry energies are smaller in high-density region 
in PM4 than in PMl, The parameter-set PM4 gives rise to a smaller effective mass than PMl, 
and the decrease of the eff'ective mass enlarges Esp-sym if the coupling strengths Ca and Ct are 
fixed. Since we fix the spin-symmetry energy at the saturation density, however, the coupling 
constant Ct becomes small for the small effective mass; then contributions from the axial-vector 
channel get larger in PM4. 

In the infinite density limit the effective mass M* goes to zero (Fig. lb) and the scalar- 
density ps converge to the finite value in the RMF theory. Thus pA and pT must have density- 
dependence similar to pB and ps, respectively; namely, when — > oo, pT converges to a finite 
value and pA is proportional to the baryon density ps- In this limit, then, only the kinetic 

^The spin-symmetry energy is written as £sp-sym ~ (p|'/6m*)(l + Go) with the nonrelativistic effective mass 
m* and the Landau parameter Go- This value is given as 26 MeV by a realistic nuclear matter calculation 
and fO - 60 MeV by the Skyrme interaction 0, |l||. 
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energy and the axial vector exchange channels contribute to the spin-symmetry energy £sp-sym- 
Since < Tiz/A >~ Xg around the spin-saturated matter, the contribution from the kinetic energy 

2/3 1/3 

is proportional to in the low density limit and to in the high-density limit. On the 
other hand the contribution from the axial- vector exchange channels is proportional to CaPb 
in the high-density limit; this behavior should be the same as the density-dependence of the 
(isospin-) symmetry energy. From this fact we can easily see that, if Ca < 0, the spin-symmetry 
energy becomes negative and the spin-polarization spontaneously occurs above a certain critical 
density. 

In Fig. 3 some quantities are shown as functions of the spin-polarization parameter Xg at 
Pb = Po (dotted line), 3po (dashed line), 5po (solid line) and 6po (chain-dotted line). We give the 
energy difference from that at the spin-saturated matter AEt/A = {Et{xs) — Et{xs = 0))/^ 
(a), M* /M (b), the total spin per nucleon < Ti^/A >= pa/ Pb and the deformation parameters 
A| and Aj. 

In Fig. 3a it can be seen that above pb > 5po the value of the spin-polarization parameter 
at the energy-minimum moves from a;^ = to a finite value, whose value becomes larger as 
baryon density increases. There is a single local energy-minimum at the fixed density. Thus 
the phase transition from normal matter to the spin-polarized one is of the second order. ^ 
To occur the first order phase-transition we need an interaction energy which is negative and 
proportional to at least < T^z/A >^. The axial-vector term in the interaction energy has the role 
to give rise to a spin-polarization, and it is exactly proportional to < Tiz/A >^. On the contrary 
we can see in Fig. 3b that the effective mass is very slightly varied with the spin-polarization 
parameter, and thereby we can conclude that the scalar interaction does not largely contribute 
the spin-polarization. 

In Fig. 3c we can see that at high density the value of < Tjz/A >= paI Pb becomes slightly 
smaller than that in the fully polarized system . From eq. ( ^T|) we can easily understand 
this behavior by considering that the small effective mass reduces the expectation value of the 
nucleon spin due to the relativistic effect {M* /E* < 1). 

In Fig. 3d it can be seen that at high density the momentum distribution gets the prolate 
deformation for spin-up nucleons and the oblate deformation for spin-down nucleons; these 
deformations enhance pA ( |3l| ) and decrease px (^). 

In Fig. 4 we show the equation of states for the fully spin-polarized nuclear matter {xg = 1), 
the density-dependence of AEt/A (a), < T^z/A > (b) and the deformation parameters for the 

^It should be interesting to compare this result with the one in quark matter 
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spin-up nucleons A(C = 1) (c), with the three spin-dependent parameter-sets (SDl— 3). As 
increasing Ca the deformation becomes larger, keeping rather large value of < S^/A >, which 
in turn implies that the phase transition more easily occurs. 

When Ca < 0, the axial- vector exchange interaction plays a role to decrease the total energy 
by enlarging < Ts^jA > with the variation of and A^q). On the other hand these variations 
enhances the kinetic and scalar exchange energies; the latter is not seen to be so effective because 
the effective mass slightly changes from the spin-symmetric value (see Fig. 3b). 

In order to examine effects of the spherical symmetry breaking, in Fig. 5, we show the spin- 
symmetry energy Esp-sym using three kinds of choices (Chi, Ch2-S, Ch2-Q). Ch2-S and -Q are 
the two versions of Ch2: the spin-vector of Ch2 (eq.(|30|)) with the spherical (Ch2-S) and the 
quadrupole-deformed momentum distribution (Ch2-Q). 

In Chi the spin-symmetry energy always monotonously increases when density becomes 
larger; the reason has been given in the previous section. In Ch2-S the qualitative behavior is 
almost the same as that for Ch2-Q though the value of £sp-sym and the critical density are always 
larger than those for Ch2-Q. When Ca < 0, the spherical symmetry breaking for the spin-vector 
makes a critical effect for the spin-polarization. In addition such effects from the choice of the 
spin- vector are enhanced by the deformation of the momentum distribution (Fig. 4c). Hence 
the axial-vector correlation between two nucleons rather easily gives rise to the ferromagnetic 
state through the spherical symmetry breaking. 

In Fig. 6 we show the density dependence of < S^/A >= pa/ Pb (a) and < /JE^/A >= pr/ Pb 
(b) at = 1 using Chi (dotted line), Ch2-S (dashed line) and Ch2-Q (solid line), and the 
quadrupole deformation parameter for spin- up nucleons A| using Ch2-Q (c). The parameter- 
sets PMl and SDl are used in this calculation. 

With increase of density, the difference among three choices becomes prominent: the total 
spin per nucleon < Ti^/A > decreases for Chi and Ch2-S, while it does not become so small for 
Ch2-Q (see Fig. 7a). This behavior can be understood from the analysis at the infinite density 
limit. In this limit, as mentioned before, the effective mass approaches to the zero value, and 
thereby < T.JA >^ x,/3 for Chi (see eq.(|3)) and < T.^/A >^ Xs/2 for Ch2-S (see eq.(|lD); 
in any choice of the spin-vector the total spin, in the relativistic framework, becomes much less 
than that in the non-relativistic one. On the other hand the choice Ch2-Q gives 



< ^^/^ >^ 2(lVt-6A,) - 2(lVe-6Ax)- (43) 



In the limit of A| ^ oo and Aj — oo, < T,z/A >— > (1 -|- Xs)/2; of course this limit makes 
infinite kinetic energy so that the finally has a moderate value. Thus the deformation of 
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the momentum distribution plays a significant role to give a large value of the total spin. The 
prolate deformation of the momentum distribution recovers the reduction of the total spin in 
the high-density region. 

Furthermore it can be seen in Fig. 6b that < > decreases for Ch2-S and Ch2-Q with 

increase of density, while its value does not becomes small. As discussed in the previous section, 
< >^ 2xs/3 for Chi and < I3T.^/A >^ for Ch2-S and Ch2-Q in the infinite density 

limit (Af* — > 0). In addition we see that the value of < l3Tiz/A > in Ch2-Q is always smaller 
than that in Ch2-S, and we know here again that the deformation of the momentum distribution 
plays a important role through the reduction of < l3Tiz/A >. 

From these results we confirm that the choice of the spin-vector is very important for the 
spin-polarization, and that we have to use Ch2 instead of Chi if Ca < 0. In the Refs. |5|, |^ the 
coupling Ca becomes positive, though the interaction is not zero-range, and then Ch2 may not 
be appropriate. Then there remain Chi and Ch3; they used Chi. 

As made before, the analysis at the infinite density limit must be useful to examine it 
qualitatively. In this limit the spin-dependent part of the energy density becomes 

esD ~ ^{CA + 'iCT)pWs [inChl] (44) 
esD ~ ^^tpWs (inChS) (45) 

where we assume that the momentum-distribution holds the spherical symmetry. To get smaller 
energy in Chi than in Ch3, 

hcA+^Cr) <^Ct (46) 
9 16 

and this equation leads to a condition that Ca < l-SSC^^, which is inconsistent with the condition 
Ct < < Ca- In addition, the deformation of the momentum distribution further reduces esn 
in Ch3. 

If Ct < < Ca, then, Ch3 must become reasonable. We can suppose that the qualitative 
behavior must be similar to that of Ch2-Q. As the density becomes larger, namely, pA converges 
to a finite value, and pT increases proportionally to the density, and then the phase-transition 
occurs at a certain critical density. In this case the momentum distribution is deformed with 
the oblate shape for the spin- up (A| < 0) and with the prolate for the spin-down (Aj > 0). From 
this consideration we can expect that the total energy should be smaller in Ch3 than that in 
Chi at any density if Ct < < Ca- 

From the above results we can see the particular role of the spherical symmetry breaking 
through the choice of the spin-vector and the distortion of the momentum distribution. They 
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becomes very important in high-density region due to decreasing of the nucleon effective mass. 

As mentioned in Sec. I the spontaneous spin-polarization is expected to occur in the high- 
density region inside neutron stars, though we have studied the magnetic properties only at the 
symmetric nuclear matter. We did because we would like to clarify the discussion by avoiding 
extending the formulation to include the isovector channel. We expect to determine the spin-spin 
interactions in the relativistic framework from some experimental information in future. However 
we may have some meaning to compare results at neutron matter with that at symmetric nuclear 
matter within the formulation given in this paper. In Fig. 7 we show the density dependence 
of the spin-symmetry energy with PMl and SD(l-3) for Ch2. The thick and thin lines indicate 
results at neutron matter and at symmetric nuclear matter, respectively. We cannot see any 
significant difference between them except that £sp-sym a^t neutron matter is a little larger than 
that at symmetric nuclear matter. 

4 Summary and Concluding remarks 

In this paper we have examined a possible mechanism of the spin-polarization of nucleons and 
discussed magnetic properties of the system. In the relativistic framework there are two kinds 
of spin-spin interaction channels, the axial-vector and tensor ones, which are reduced to the 
same interaction channels in the nonrclativistic framework. If the interaction energies from two 
channels have opposite signs, there is a second-order phase transition to a spin-polarized state. 
Though the effects from two channels are counterbalanced with each other around the normal 
density, the channel with the negative sign becomes dominant, suppresses the spin-symmetry 
energy with increase of density and induces a phase-transition in a certain critical density pc- In 
this mechanism the spherical symmetry breaking through the spin vector and the momentum 
distribution plays a significant role; the spherical symmetric calculation (Chi) cannot describe 
such a phase-transition sufficiently. 

These qualitative findings can be easily derived from the analysis at the infinite density limit, 
which is equivalent to the ultra relativistic limit in the present framework. Actual numerical cal- 
culations confirm them, so that the consideration in this limit is very useful to predict qualitative 
behaviors of matter at high density. 

If the tensor channel is largely negative, the total nucleon spin becomes large while the 
total spin converges to zero at the infinite density limit. Thus such system cannot make strong 
ferromagnetism even at high density. Inside actual neutron stars the density is not infinite, and 
we has not known how large ferromagnetism is necessary to explain magnetars. Then we cannot 
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deny this possibility though we omit it at the present calculations. We should keep it in mind 
for future. 

In this work we represent the strength of these channels with the couplings Ca and Ct by 
using the zero-range approximation, which makes the Fock exchange interaction local. Then 
we find that, if Ca < —100, the phase-transition occurs in the reasonable density pB < 5pQ, 
which can be realized inside neutron stars. However, the main purpose of this paper is to 
reveal the characteristic features of the relativistic ferromagnetism of nuclear matter within the 
RHF framework. If we want to get a realistic conclusion about the critical density, we need 
to determine two couplings Ca and Ct individually from the experimental information, while 
experiments can at present give only nuclear spin properties in low density region around and/or 
below the normal density. Furthermore, we need the isovector interactions when we consider 
neutron-star matter. 

Two approximations have been introduced in this work. One is that we have applied a 
variational approach, avoiding the complete self-consistent calculation. If we can do it, some 
points are improved; the single particle energies become different between the spin-up and - 
down states, which should in turn determine the modification of the Fermi sea. Because of 
the variational principle, anyway, these improvements must reduce the total energy in the spin- 
polarized system. Then the phase transition may occur in the lower density and/or with the 
smaller couplings of Ca than those in the present calculations. 

The other is that we have used the zero-range approximation for the meson propagators and 
discarded the finite-range effects in the nucleon interactions; the zero-range approximation is 
equivalent to neglect the momentum-dependence of the self-energies and should be reasonable, 
at least, at low densities s.t. nia ^ pf- As for the spin-independent parts it has been reported 
that the momentum-dependence of the self-energy has a role to reduce the total energy per 



nucleon [Ex /A) in the high-density region [^, 16| and to suppress largely the Fermi velocity 



|17| , particularly in the low density region. However these effects does not affect the nuclear 
equation of state at zero temperature. As for the spin-dependent parts the situation must be 
similar. 

These approximations used in our calculations must never spoil our qualitative findings of 
the relativistic ferromagnetism. 

We do not clearly know what phase actually appears inside neutron stars. A quark matter 
1^] is one possibility and the hyperonic matter is also possible. In the latter case we need to take 
into account the interaction between hyperon and nucleon and that between hyperons. 
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We compare here our results with the previous one given by one of the authors (T.T.) for 
quark matter |^]. It has been shown that ferromagnetism of quark matter may occur at low- 
density region and it should be the first-order phase transition; all the quark spins suddenly align 
at the critical density and there is no partially polarized state. The most important difference 
from nuclear matter is that quarks interact with each other through only the vector interaction 
by gluons and there is no direct interaction due to color neutrality of quark matter. Moreover 
there is no tensor channel and the gluon propagator corresponds to the zero-mass limit of A^^. 
As a result "Dyiiq) given in Eq. (5) is always negative and I't(q) = 0, which is most favorite 
situation for the system to be ferromagnetic. In this case the momentum dependence of the 
propagator becomes essential and the phase transition is of the first order due to this effect; if 
we replace the propagator by some constant, we shall see that phase transition becomes of the 
second order and there is a partially polarized state as is seen in this paper. 

Our findings of the momentum dependence of the spin orientation and possible deformation 
of the Fermi sea have not been taken into account in ref.Q. It may be interesting to see these 
effects in the context of ferromagnetism of quark liquid |18]. 
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Figure Captions 



Fig. 1 Density-dependence of the total energy per nucleon (a), and the ratio of the effective 
masses to the bare masses for nucleon (b). The solid and dashed lines indicate the results 
for PMl and PM4, respectively. 

Fig. 2 Density-dependence of the spin-symmetry energies {ssp-sym) with PMl (a) and PM4 (b). 
The long-dashed, dashed and solid lines indicate results for Ca = (SDl), = —50 (SD2) 
and = —100 (SD3), respectively. 

Fig. 3 Energy difference between the spin-polarized and saturated systems (a) , the effective mass 
M* normalize by the bare nucleon mass M (b) the total spin per nucleon < T>z/A > (c). 
and the deformation parameter A (d) versus the spin-polarization parameter {p-[ — Pi)/ Pb- 
The dotted, dashed, solid and chain-dotted line indicate results at pB = Po, Pb = ^Po, 
Pb = 5po and Pb = Qpo, respectively In the third column (c) the thick and thin lines 
indicate ones for the spin-up and spin-down, respectively. 

Fig. 4 The density-dependence of AEt/A (a), < Sz/A > (b) and the deformation parameters 
for the spin-up nucleons \{( = 1) (c) with SDl (dotted line), SD2 (dashed line) and SD3 
(solid line), in the fully spin-polarized nuclear matter {xg = 1). 

Fig. 5 Density-dependence of the spin-symmetry energies {sspsym) with PMl calculated with 
the method Chi (a), Ch2-S (b) and Ch2-Q (c) (seeing text). The meaning of the lines are 
shown in Fig. 2. 

Fig. 6 The density-dependence of < az/A > (a) < Paz/ A > (b) and the deformation parameters 
for the spin-up nucleons \{( = 1) with SDl (dotted line), SD2 (dashed line) and SD3 
(solid line) in the fully spin-polarized nuclear matter (xs = 1)- For the spin- independent 
parts the parameter-sets of PMl are used. 

Fig. 7 Density-dependence of the spin-symmetry energies {ssp-sym) with PMl. The long-dashed, 
dashed and solid lines indicate results for Ca = (SDl), = -50 (SD2) and = -100 (SD3), 
respectively. The thick and thin lines indicate results at the neutron matter and at nuclear 
matter, respectively. 
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5<T gui B„ 


PMl 


9.408 9.993 23.52 5.651 


PM4 


11.05 12.64 18.89 7.158 



Table 1: Parameter sets for the RH calculation in this paper. In all cases have used 777.^= 550 
MeV, m^= 783 MeV and = 0. 





PMl 


PM4 




Ca Ct 


Ca Ct 


SDl 


9.993 


1.1365 


SD2 


-50 56.10 


-50 48.32 


SD3 


-100 104.5 


-50 96.27 



Table 2: Parameter sets for the spin-dependent interactions in this paper. 
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